Newton’s method’s basin of attraction for sign-changing solutions of concave and convex nonlinearities

In this paper, we will illustrate and visualize the chaotic behavior for solving the sign-changing solutions of concave-convex nonlinearities. We propose contractor iteration algorithm and bifurcation method to improve the computation speed and avoid local convergence of Newton-Raphson method. By this method, sign-changing numerical solutions will be computed and visualized in various domains.     

An improved method for shape measurement using two-dimensional digital image correlation

Shape measurement is a significant application of digital image correlation (DIC). An improved method that combines a rotatable plane mirror is proposed to measure the shape of an immovable object. In this method, two images, one before and the other after rotating the plane mirror, are obtained and then in-plane translation which related to the shape of the detected object can be calculated by the use of two-dimensional digital image correlation (2D DIC). The relationship between the in-plane translation and the shape of the object is described. Experimental results show that the proposed method is feasible for shape and distance measurement with high accuracy.     

A Combined Relaxation Method for a Class of Nonlinear Variational Inequalities

In this paper, we describe a class of combined relaxation methods for the non strictly monotone nonlinear variational inequality problem, which involves a max-type convex function. This method is readily implementable and attains a linear rate of convergence under certain additional assumptions.     

Monte Carlo EM加速算法

EM算法是近年来常用的求后验众数的估计的一种数据增广算法, 但由于求出其E步中积分的显示表达式有时很困难, 甚至不可能, 限制了其应用的广泛性. 而Monte Carlo EM算法很好地解决了这个问题, 将EM算法中E步的积分用Monte Carlo模拟来有效实现, 使其适用性大大增强. 但无论是EM算法, 还是Monte Carlo EM算法, 其收敛速度都是线性的, 被缺损信息的倒数所控制, 当缺损数据的比例很高时, 收敛速度就非常缓慢. 而Newton-Raphson算法在后验众数的附近具有二次收敛速率. 本文提出Monte Carlo EM加速算法, 将Monte Carlo EM算法与Newton-Raphson算法结合, 既使得EM算法中的E步用Monte Carlo模拟得以实现, 又证明了该算法在后验众数附近具有二次收敛速度. 从而使其保留了Monte Carlo EM算法的优点, 并改进了Monte Carlo EM算法的收敛速度. 本文通过数值例子, 将Monte Carlo EM加速算法的结果与EM算法、Monte Carlo EM算法的结果进行比较, 进一步说明了Monte Carlo EM加速算法的优良性.     

掠入射较高功率密度驱动Ni-like AgX射线激光的理论研究

 用系列程序对掠入射方式驱动亚稳态类镍银碰撞激发机制进行了研究。继低功率密度驱动获得有效增益系数为23.8 cm^-1结果之后，在4.3×10¹³ W·cm^-2较高功率密度驱动时，增益区移到驱动激光折返点附近，有效增益系数达到40 cm^-1以上，比正入射提高了2.4倍。用50 J激光能量驱动单靶，就能达到有效增益长度积为40的深度饱和增益。     

关于剩余类环Zn和Galois环GR（p^e，n）上的理想双同态门限体制

本文讨论了两类特殊环上的理想双同志门限体制．利用门限体制的完善性及同态性，证明了以剩余类环Zn，和Galois环GR(p^e，n)为密钥空间的理想双同志(t，ι)-门限体制是平凡的．     

Global analysis of continuous analogues of the Levenberg-Marquardt and Newton-Raphson methods for solving nonlinear equations

Summary Global analyses are given to continuous analogues of the Levenberg-Marquardt methoddx/dt=−(J ^t(x)J(x)+δI)⁻¹J^t(x)g(x), and the Newton-Raphson-Ben-Israel methoddx/dt=−J ⁺(x)g(x), for solving an over- and under-determined systemg(x)=0 of nonlinear equations. The characteristics of both methods are compared. Erros in some literature which dealt with related continous analogue methods are pointed out. The Institute of Statistical Mathematics     

Pressure and time treatment for Chebyshev spectral solution of a Stokes problem

To investigate the influences of time scheme, pressure treatment and initial conditions in incompressible fluid dynamics, a Stokes problem is solved numerically on a slab geometry within the framework of spectral approximation in space. Four algorithms are examined: splitting schemes, influence matrix method, penalty formulation and pseudo-spectral space-time technique. It is shown that splitting schemes are less accurate than the other processes. Furthermore, the initial field should respect a compatibility condition to avoid singularities at the initial time. If it is not possible to build such a compatible field, the numerical procedure has to present good damping properties at the first steps of the time integration.     

A finite element convergence study for shear-thinning flow problems

The solution of the non-linear set of equations arising from the application of the finite element method to non-Newtonian fluid flow problems often requires large amounts of computer time. Four iteration schemes (Picard, Newton-Raphson, Broyden and Dominant Eigenvalue method) are compared in three different flow geometries using a shear-thinning fluid model. Points of comparison involve the computer time necessary to converge the equations, ease of implementation, radius of convergence and rate of convergence.     

Solution of three-dimensional natural convection about a vertical square rod in a cylindrical enclosure

The problem of three-dimensional laminar natural convection in a vertical enclosure with an inner square rod is treated by a numerical method in boundary-fitted co-ordinates. The inner and outer cylinders are heated and cooled, respectively, to maintain different constant surface temperatures. The horizontal enclosure surfaces are maintained at adiabatic conditions. The Prandtl number is that of air, 0·703, and the Rayleigh numbers span the conduction, transition and boundary layer regimes of flow. The radius ratio is 1, and the aspect ratio (cylinder length divided by maximum annular gap) is 1. The results of the study provide data useful in the design and performance assessment of nuclear reactor spent fuel shipping casks.